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Abstract. Current cosmological models and data suggest the existence of a Cold Dark Matter (DM) 
component, however the nature of DM particles remains unknown. A favored candidate for DM is 
a Weakly Interacting Massive Particle (WIMP) in the mass range from 50 GeV to greater than 10 
TeV. Nearby dwarf spheroidal galaxies (dSph) are expected to contain a high density of Dark Matter 
with a low gamma-ray background, and are thus promising targets for the detection of secondary 
gamma rays at very high energies ( VHE, E > 0. 1 TeV) through the annihilation of WIMPS into SM 
particles. Presented here are recent VERITAS observations of dSph, including a deep exposure on 
Segue 1. Limits are derived for various annihilating and decaying dark matter particle models. 
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INTRODUCTION 

The indirect search for very high energy (VHE) y-rays resulting from the annihilation 
of Dark Matter (DM) particles into SM particles provides an important complement 
to that of direct searches for DM interactions and accelerator production experiments. 
Among theoretical candidates for the DM particle [1], a weakly interacting massive 
particle (WIMP) is well-motivated since it naturally provides the measured present day 
cold DM density [2]. Candidates for WIMP dark matter are present in many extensions 
of the SM of particle physics, such as supersymmetry (SUSY) [3] or theories with extra 
dimensions [4]. In such models, the WIMPs either decay or annihilate into standard 
model particles, producing either a continuum of y-rays with energies up to the DM 
particle mass, or monoenergetic y-ray lines. 

Promising targets for indirect DM searches are nearby regions with high DM density. 
The Galactic Center is likely the brightest source of y-rays resulting from DM annihi- 
lations, however the detected VHE y-ray emission is coincident with the supermassive 
black hole Sgr A* and a nearby pulsar wind nebula [5], motivating searches for DM 
annihilation in the Galactic Center halo where the VHE y-ray background is expected 
to be significantly lower [6]. Alternatively, dwarf spheroidal galaxies (dSphs) of the 
Milky Way are attractive targets for DM searches. This is due to their relatively well- 
constrained DM profiles derived from stellar kinematics, proximity of 20-100 kpc, and 
general lack of active or recent star formation, suggesting a relatively low background 



http://veritas.sao.arizona. edu/ 



TABLE 1. VERITAS observations of dwarf spheroidal galaxies (dSph). 



dSph 


Dist. 


Obs. Period 


Exp. 


Sig. 


UL 


E Th 


Flux UL 




(kpc) 




(hr) 




(cts) 


(GeV) 


(cm 2 s _1 ) 


Draco 


80 


2007 Apr-May 


18.4 


-1.5 


18.8 


340 


0.49 


Ursa Minor 


66 


2007 Feb-May 


18.9 


-1.8 


15.6 


380 


0.40 


Willman 1 


38 


2007-2008 


13.7 


-0.1 


36.7 


320 


1.17 


Bootes 1 


62 


2009 Apr-May 


14.3 


1.4 


72.0 


300 


2.19 


Segue 1 


23 


2010-2011 


47.8 


1.4 


102.5 


300 


0.80 



from conventional astrophysical VHE processes [7]. Here we present results for a search 
of DM annihilation signal from five dSphs observed with VERITAS. 

OBSERVATIONS 

VERITAS is an array of four imaging atmospheric Cherenkov telescopes (IACTs), 
each 12 m in diameter, located at the Fred Lawrence Whipple Observatory in southern 
Arizona, USA. Each VERITAS camera contains 499 pixels (0.15° diameter) and has a 
field of view of 3.5°. In the summer of 2009 the first telescope was moved to its current 
location in the array to provide a more uniform distance between telescopes, improving 
the sensitivity of the system [8]. VERITAS is sensitive over an energy range of 100 
GeV to 30 TeV with an energy resolution of 15%-20% and an angular resolution (68% 
containment) of less than 0.1° per event. A source with a 1% Crab Nebula flux can be 
detected by VERITAS in approximately 25 hours. 

Table 1 lists the VERITAS observations of dSphs. Data-quality selection requires 
clear atmospheric conditions, based on infrared sky temperature measurements, and 
nominal hardware operation. The total exposures after dead time correction for each 
dSph are listed in table 1 . The observations of Draco, Ursa Minor, Willman 1, and Bootes 
1 were taken before the first telescope was moved, while the Segue 1 data was taken in 
the current array layout. All data were taken during moon-less periods in wobble mode 
with pointings of 0.5° offset [9]. Data reduction followed the standard methods, yelding 
consistent results with two analysis packages [10]. 

FLUX UPPER LIMITS 

As listed in table 1, no significant VHE y-ray signal is detected from the dSph observed 
by VERITAS. Upper Limits (UL) on the integrated flux are calculated based on the 
assumption that the dominant source of y-rays is DM annihilation. The y-ray differential 
energy spectrum from DM particle annihilation depends on the DM model, in particular 
on the branching ratios to the final state particles. In almost every channel (excepting the 
and channels), the y-ray emission mostly originates from the hadronization 

of the final state particles, with the subsequent production and decay of neutral pions. 
Annihilation to three different final state products is considered independently of the 
dark matter model, in each case with a 100% branching ratio. Figure 1 (left) shows 
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FIGURE 1. Left: Dark Matter annihilation spectra for four different final state products (W + W~, Z°Z°, 
bb, and t + t~), extracted from PYTHIA 8.1. Right: Upper limits (95% CL) on the integrated y-ray flux 
above E m i„ = 300 GeV from the VERITAS observations of Segue 1 considering dark matter particle 
annihilation/decay for three different channels: bb, and T + T~. [11]. 



the annihilation spectra for the channels bb, and t + t~. For each channel, the 

y-ray spectrum has been simulated with the particle physics event generator PYTHIA 
8.1. Figure 1 (right) shows the ULs (95% CL) on the integrated flux above E m i n = 300 
GeV from VERITAS observations of Segue 1 as a function of the DM particle mass. 
The integrated flux ULs range between 0.3% and 0.7% of the Crab Nebula integral flux, 
depending on the DM particle mass. 



UPPER LIMITS ON THE DM ANNIHILATION CROSS-SECTION 

The y-ray flux from the annihilation of DM particles in a spherical DM halo is deter- 
mined in part by J(Afl), the squared DM density integrated along the line of sight and 
over a solid angle given by the size of the signal search region defined in the VERI- 
TAS analysis. The estimate of J(A^) requires a modeling of the dSph DM profile, which 
for the observations of Segue 1 an Einasto profile [12] is applied, resulting in J(A^) 
= /. / x 10 18 GeV 2 cm 5 sr. The upper limit (UL) on the total velocity-weighted annihi- 
lation cross-section (c7v) 95%CL is then calculated using the UL on y-like events Ny cCL , 
the total observation time T j, s and the effective area A e ff for the dataset, and the DM 
mass m^: 

<*v) 95%CL = x T CL 1. (1) 



The main uncertainties in the DM limits computed here are from the modeling of the 
DM density profile. Unlike the classical dSphs, there is no high-statistics star sample 
for Segue 1 . This prevents an accurate modeling of the DM distribution. Assuming an 
Einasto profile [12], the systematic uncertainties in the dark matter profile modeling can 
change the astrophysical factor, and hence they scale up or down the limits by a factor 
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FIGURE 2. Upper limits (95% CL) from VERITAS observations of Segue 1 on the WIMP velocity- 
weighted annihilation cross-section (<7v) as a function of the WIMP mass, considering different final state 
particles. The grey band area represents a range of generic values for the annihilation cross-section in the 
case of thermally produced dark matter. Left: hadronic channels bb, and t + t~. Right: leptonic 

channels e + e~ and jX + jJ.~ [11]. 



of 4 at the 1 o level [13]. Future spectroscopic surveys should increase the Segue 1 star 
sample and eventually reduce the uncertainties on its dark matter content. 
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